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Abstract 

We calculate the three-point correlation function of the comoving curvature pertur- 
bation generated during an inflationary epoch driven by false vacuum energy. We get a 
novel false vacuum shape bispectrum, which peaks in the equilateral limit. Using this 
result, we propose a scenario which we call "old curvaton". The shape of the resulting 
bispectrum lies between the local and the false vacuum shapes. In addition we have a 
large running of the spectral index. 
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1 Introduction 



Inflation [TJ [2] nowadays is accepted as the standard paradigm to solve many cosmological prob- 
lems, and to provide the appropriate initial conditions for the hot big bang cosmology. One of 
the strongest supports of this early era of inflation is the observed primordial curvature pertur- 
bation 1Z C , whose power spectrum V-r is nearly scale invariant with its index being n-n pa 0.96, 
with almost perfect Gaussian statistics [3]. Conversely, these observations have been acting 
as a powerful discriminator among the models of inflation based on particle physics [I] and 
placing strong constraints on the realization of inflation, e.g. in the context of super gravity [5]. 
In the next few years, we will obtain more precise data on the primordial density perturba- 
tions by ongoing and forthcoming experiments like the Wilkinson Microwave Anisotropy Probe 
(WMAP) and the Planck satellites, and thus will be able to constrain the models of inflation 
even stronger than now. 

In the light of these upcoming observations, there have been increasing interests in the 
non-Gaussian signature of the primordial perturbations. Currently, the non-Gaussianity of the 
cosmic microwave background (CMB) fluctuations is constrained to be |/nl| ^ 100 [3], with 
/nl being the non- linear parameter [6], which translates into less than 0.1% of non- Gaussian 
contribution. The sensitivity is expected to be improved at the level of /nl ~ on the CMB 
scales [TJ, and indeed primordial non-Gaussianity will be another powerful probe of inflation: 
if substantial non-Gaussianity were to be detected, we should look for new models beyond the 
simplest single field slow-roll inflation models, where /nl 0(1) jS]. A short list of the models 
where large non-Gaussianity may arise includes Dirac-Born-Infeld (DBI) inflation [9], single 
field inflation with featureaj [H], multi- field inflation^ [13], curvaton [H] , and so on. 

In this note, we study a new possibility of generating large amount of highly non-Gaussian 
curvature perturbation: we consider a period of inflation supported by a non-zero vacuum 
energy. This is in fact the original model of inflation [TJ and has received a renewed interest 
in the context of vast string landscape [15]. But 1Z C generated during this stage by vacuum 
fluctuations has been unknown until very recently [16]. The difficulty is that while the inflaton 
field is well anchored in a false vacuum, the comoving curvature perturbation, which is given 
by 



is ill defined. The origin of this difficulty is that during false vacuum inflation the inflaton <j) is 
classically not moving, and there is no preferred time direction to define comoving hypersurfaces. 
But this phase should be broken, and indeed in Ref. [16] Vn was calculated by incorporating a 
regularizing mass which scales as a~ 2 . The results are such that V-r. is dependent on the scale 
of the regularizing mass, with steeply blue index n-ji = 4 [16] and nearly x 2 distribution^) |17j . 
Thus, we need an additional mechanism to generate 1Z C to match observations. 

In addition to false vacuum inflation, we incorporate the curvaton mechanism [19] to gener- 

4 In this case, we also have characteristic modulation and oscillation in the power spectrum [ID] . 
5 Note also an intermediate model of quasi-single field inflation [12] . 

6 Note that in Ref. 18. a similar situation was discussed, with the power spectrum induced from the false 
vacuum inflation and the quantum tunneling in spatially open universe. We expect similar results, i.e. Vn is 
suppressed on large scales and is highly scale dependent with n-ji = 4, since these are related to the properties 
of heavy fields in de Sitter space." 
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ate the Gaussian contribution of TZ C . In this combined scenario of false vacuum inflation and the 
curvaton mechanism, or shortly "old curvaton" scenario, as we will see we can build a concrete 
model with observationally consistent Vn, while providing substantial degree of non-Gaussian 
signature. Especially, the bispectrum lies between the so-called equilateral and local types, 
giving a concrete example of mixed non-Gaussianity. The structure of this note is as follows. 
In Section [2j we explicitly calculate the three-point correlation function, or the bispectrum of 
1Z C using the results of Ref . [16] . Then in Section [3] we incorporate the curvaton mechanism 
and present the resulting power spectrum and the bispectrum in the old curvaton scenario. We 
conclude in Section HI Some technical detail is given in Appendix. 



2 Bispectrum from false vacuum inflation 

We consider a theory where the inflaton sector is described by the Einstein gravity with scalar 
Lagrangian, 

77? 2 I 

£ = -fR - 2<W# - v{<t>) , (2) 

where mpi = (SttG)^ 1 ^ 2 ~ 2.4 x 10 18 GeV. The potential is assumed to have the form 



V(<P) = V + 1 -(-ml + ^)<p 2 . (3) 



We are interested in the epoch while m 2 s = — m 2 + /J 2 /a 2 > so that cf) is well anchored 
at <p = 0. During this stage, as shown in Appendix |A], we can find the comoving curvature 
perturbation TZ C as 

1 f d 3 a 

K < « = 2a^ + P ) re „( W ) 2 J [ " + q ' (<1 - k)1 ' <4) 

The convolution in the momentum space will lead to a loop integral if we calculate the two- 
point or three-point correlation function. It's worth to note that this kind of loop integral is 
not the real loop integral in interactional quantum field theory. It is the so-called c-loop [20] ■ 
As can be read from (TJJ, the curvature perturbation is quadratic in Scf), so we need not follow 
the conventional approach [8]. The bispectrum arise from c-loop integration between different 
modes of curvature perturbation. The three-point correlation function of TZ C is explicitly written 
as 

-qi ^q2 0k 2 -q2 0q3 0k 3 -q 3 ) 



(5) 



[2a 2 (p + p)Uw) 2 } 3 J ( 2ff ) 3 ( 27r ) 3 ( 27r ) 3 
x [V 2 + Qi ■ (qi - ki)] [/i 2 + q 2 ■ (q 2 - k 2 )] [fi 2 + q 3 • (q 3 - k 3 



We can after some straightforward calculations find the expectation value by explicitly plugging 
(1371) into but we can anticipate the result by considering the possible contractions between 
0's. It is obvious that there are 8 possible physically meaningful contractions: TZ c (ki) to 7£ c (k 2 ) 
or lZ c (k 3 ), so we have 2 choices. Further, since each 7£ c (kj) has 2 field contents, there are 2 
choices for each: one of the fields in (say) TZ c (ki) can be correlated to one in (say) 7£ c (k 2 ) and 
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the other in 7£ c (ki) can be correlated to one in TZ c (k 3 ). This is to avoid self-contraction of 
lZ c (kz). Thus total 2x2x2 = 8 different ways of contraction. These are explicitly written as 

i — I i i "I 1 i — I n i n I 

((0k 1 0k 1 -q 1 )(0k 2 0k 2 -q 2 )((/'k3 0k3-q3)) + ( (0k a 0ki -qi ) (</>k 2 </>k 2 -q 2 ) (</>k 3 </>k 3 -q 3 ) ) 



i i i i i i i i i i n 

+ ((0k 1 0k 1 -q 1 )(0k 2 0k 2 -q 2 )(0k3 0k3-q 3 )) + ((^k^-qj (</>k 2 </>k 2 -q 2 ) (</>k 3 </>k 3 -q 3 ) ) 

i i i r= > i i i I i i I 

+ ((</>k 1 0k 1 - qi )(0k 2 0k 2 -q 2 )(^k3 0k3-q 3 )) + ( (<£k a 0k x - qi ) (0k 2 0k 2 -q 2 ) (</>k 3 </>k 3 -q 3 ) ) 



I II I I I I l l 



+ ((0k 1 0k 1 -q 1 )(0k 2 0k 2 -q 2 )(0k30k3-q 3 )) + ( (0k x 0k x -q J (0k 2 0k 2 -q 2 ) (0k 3 0k 3 -q 3 ) ) , (6) 

and these terms exactly correspond to the terms with (^(ki + k2 + k 3 ) which we can obtain 
by direct computations. Indeed, after some trivial calculations, we find that ([5]) can be written 

as 

(^ c (k 1 )^ c (k 2 )^ c (k 3 )) 
4^ 3 )(k! + k 2 + k 3 ) 

x J ^{I^H^I^-qlH^+qll 2 [/i 2 + q- (q-kx)] [fi 2 + q • (k 2 + q)] [/i 2 + (q - k x ) ■ (k 2 + q) 

+ |^| 2 |^|k 1 -q|| 2 |^|k3+q|| 2 + q ■ (q - ki)] [ll 2 + q • (k 3 + q)] [ll 2 + (q - kx) ■ (kg + q)] } 

(7) 

Where ip k is the mode function of perturbation 0k (see (A. 11) for example), 

0k = ctk^fc + a>- k ¥l ■ ( 8 ) 

Since the two terms of the integrand of (J7J) are exactly the same with the permutation of fc 2 
and &3, we only consider the first term. 

As can be read, is a very complex function of three momenta. However, since we are 
interested in large scales, we can expand each term of (jTJ) and only keep the lowest non-zero 
contribution. Then, explicitly using the asymptotic form of the mode function (p k when \rj\ ^> 1, 

2 — ► \ n 9 
^ ' H»i 2^/WTp? 

and integrating over angles, up to linear order of k we can find 
J rffi|^| 2 |^| kl _ q || 2 |v?| k2+q || 2 [ll 2 + q • (q - ki)] [ll 2 + q • (k 2 + q)] [ll 2 + (q - kj) ■ (k 2 + q 

= l(Hv)*tf + S) 3/ \ (10) 

and k dependence appears only beyond quadratic order. Now we have to perform the integration 
with respect to the magnitude q — |q|. If one naively integrates from to infinity, the integral 
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badly diverges. However, actually we should appropriately regularize the integral as follows: 
we are interested in the large scale curvature perturbation produced during the false vacuum 
inflation stage. Thus, we can trace the momentum up to the scale at which false vacuum 
inflation ends. That is, we should consider up to the scale fc* which crosses the horizon at the 
moment the potential becomes concave down. This can be found from the potential fl3]): it 
becomes concave down when the effective mass is zero, i.e. m 2 = fi 2 /al, which in turn gives 
the scale factor at this moment as a* = fi/m^. The momentum scale which crosses the horizon 
at this moment is thus = a*H = (H/m^/j,. The modes with larger wavenumbers k > k* exit 
the horizon after the false vacuum inflation phase. So for the purpose of bispectrum produced 
during false vacuum inflation, one should cutoff the integration at k = k*. 
After imposing the cutoff at k*, the integration becomes 

J rf 3 g|^| 2 |^| kl _ q | | 2 |^| k2+q | | 2 [fJ 2 + q • (q - k^] [/J 2 + q • (k 2 + q)] [fx 2 + (q - k x ) ■ (k 2 + q 



7T 



96 
where 



fi- 
ll 



(11) 



k 2 ( k 2 kt 

l + 3 + 14-^ + 8^ 



3 log 




(12) 



For example, when = H, we have f(k*/fi) m 30.6318. Thus, with the other term of (J7J) for 
which we can simply replace £; 2 by &3, the full three-point correlation function of 1Z C is found 
to be, up to an 0(1) constant, 



(^ c (k 1 )^ c (k 2 )^ c (k 3 )) = 5® fa + k 2 + k 3 )^-|^/ (- 

6A m^fi v fi 

where we have used (|3"T|) . Comparing this with the definition of the bispectrum 
(^ c (k 1 )^ c (k 2 )^ c (k 3 )) = (2n) 3 S^(k 1 + k 2 + k 3 )^(k 1 , k 2 , k 3 ) , 



(13) 



(14) 



we can find the bispectrum of the curvature perturbation 1Z C produced during false vacuum 
inflation as 



B w (ki,k 2 ,k 3 



(2tt) 4 H 6 (K 



3A 3 mS" 6 ' 



(15) 



Note that B-R,(ki, k 2 , k 3 ) is independent of kj at linear order of ki/fi. This is very different from 
non-Gaussianity produced by any other known models. 

In Fig. [H we show the dimensionless shape function (/c 1 fc 2 fc 3 ) 2 i?7^(k 1 , k 2 , k 3 ) as a function of 
kijk\ and k^jk\. It is clear that the bispectrum exhibits its maximum value at the equilateral 
limit k\ = k 2 = k 3 . This means, unlike the so-called local type non-Gaussianity where non- 
linearity arises due to the classical super-horizon evolution of the curvature perturbation, the 
curvature perturbation itself is intrinsically highly non-Gaussian. Indeed, as we can see from 
(jl]), 1Z C oc 2 with being nearly Gaussian, thus TZ C follows x 2 statistics. A more close study 
will follow separately [T7] . 
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Figure 1: The shape of the bispectrum produced during false vacuum inflation. We have the 
maximum value of the bispectrum when all the momenta are of comparable magnitude. 

3 Old curvaton scenario 

In the previous section, we calculated the bispectrum of 1Z C produced during false vacuum in- 
flation. However, false vacuum inflation cannot reproduce the observed nearly scale invariant 
curvature perturbation with nearly Gaussian distribution. To generate the Gaussian contribu- 
tion, we introduce a curvaton field o. As we will see shortly, in this "old curvaton" scenario 
we can obtain the desired nearly scale invariant power spectrum of the Gaussian curvature 
perturbation, along with novel shape of the three-point correlation function as well as large 
running of the spectral index. 

We employ the standard curvaton scenario for producing the nearly scale invariant power 
spectrum. The curvaton a has very small mass, no direct coupling with the inflaton field, and 
sub-dominate energy density during inflation. After inflation, the curvaton field eventually 
dominates the energy density either during curvaton oscillation stage or a secondary inflation 
stage. Then the isocurvature perturbation stored in the curvaton field is converted to the 
curvature perturbation. 

In the following calculation, we will denote as IZ^ and lZ a the curvature perturbations on 
the comoving slices of the inflaton (including possible radiation component) and the curvaton, 
respectively, and use 1Z C to denote the curvature perturbation on the comoving slice of the 
total energy- momentum tensor. As shown in Ref. [2TJ, TZ^ and lZ a are separately conserved 
quantities as long as there is no direct coupling between the inflaton and the curvaton er, 
which is assumed in our scenario as well as in the standard curvaton scenario. 

The comoving curvature perturbation defined from the total energy-momentum tensor can 
be written as 




(16) 
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where the constant r is given by 

Po 



(17) 

dec 



P<j> + Pa 

which is evaluated when curvaton decays: after then, TZ C becomes constant. From ffl6l) . the 
power spectrum can be written as 

Vn = {l-rfVn,+r 2 Vn a , (18) 

where we have assumed no cross correlation between and a. Actually, such correlation can 
be induced by the gravitational coupling. But this cross correlation can be absorbed into the 
redefinition of IZ^ and TZ a without modifying any other calculation present in this paper: we 
will present more discussion on this redefinition at the end of this section. The spectral index 
is 

a log k V n 
where n a = dlogVn a /dlogk. The running of the spectral index is written as 

a n = = {n n - 1) - {n n - l) 2 + 6 + [a a - {n a - 1) + {n a - l) 2 - 6] r 2 ^ , (20) 

a log A; P-R 

where a a = dn a /dlogk. Note that the running of the index is of the same order of magnitude 
as the spectral index. We can obtain a large running in the old curvaton scenario. 

The total three-point correlation function in the old curvaton scenario can be calculated as 

(^ c (k 1 )K(k 2 )^ c (k 3 )) = (1 - r) 3 (^(k 1 )^(k 2 )^(k 3 )) + r 3 (^(k 1 )^ (T (k 2 )^ CT (k 3 )) , (21) 

with (TZ^(ki)7l c j ) ('k2)7l ( f > (k 3 )) given by (TT3]) . (7?. (J (ki)7?. cr (k2)7?. (J (k3)) has a typical local shape, 
which is standard in the curvaton scenario written as[14j 

r 3 (^(k 1 )^(k 2 )^(k 3 )) = (2vr) 7 5 3 (k 1 + k 2 + k 3 ) (-^K) ~ | " f ) ^ • (22) 

Depending on which term is dominant in (12 the total three-point correlation function in 
the old curvaton scenario can exhibit different shape: if the inflaton contribution is larger, the 
shape coincides with strong equilateral type shown in Fig. [TJ while the curvaton is dominating 
we have the well-known local shape as shown in Fig. |2J More generally, in the old curvaton 
scenario in general we have a mixed shape, lying between these two extremes. An example of 
mixed shape is illustrated in Fig. |3j 

We can also match our mixed shape non-Gaussianity with the local shape at k\ = k% = k 3 
to obtain a non-Gaussian estimator as 

/NL "4r 3 6 1 } 27A*V%mltjfi f \») " 1 ' 

Using the expression of Vn^ [IS], ( 123]) can be further simplified to 

fmix _A_5_5r r? lOA^> (K\ (VnA 2 
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Figure 2: The local shape of bispectrum produced by the nonlinearity of curvaton. The bis- 
pectrum is maximum when either A; 2 or /c 3 is comparable to k\, while the other one negligible. 




Figure 3: An example of the mixed shape between the local and the false vacuum shapes. The 
momentum dependence of the bispectrum is of the form Bn(ki, k 2 , k 3 ) ~ 20 + X^^/TL^; 
where the former and the latter are the false vacuum and the curvaton contributions, respec- 
tively. 
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where B = 0(0.1), and m$ is a free parameter in our scenario. We can get a positive and large 
/nl by tuning m$ 3> i? 2 . Note that /^l* is strongly scale dependent: the scale dependence is 

encoded in the fraction (Vk^/V-r,) 2 ■ 

Finally, we would like to make some remarks on the aforementioned redefinition of the cur- 
vature perturbations. In the calculation, we have ignored the gravitational coupling between 
4> and a, and assume these two fields are uncorrelated. Actually, they should interact gravita- 
tionally so that part of 5a is induced from via 71$ at second order. However, the curvature 
perturbation from this part of 5a is proportional to 71$, thus can be absorbed into the redefined 
7Z$. To be explicit, one can decompose 7Z a into two parts, 

n a = n { j ac) + n£ Qduced) , (25) 

where TZ^^ originates from the vacuum fluctuation of the curvaton, thus has nearly Gaussian 
distribution, while 7^^ induced ) j s induced by the gravitation potential during the false vacuum 
stage of inflation, thus follows x 2 distribution. Then the total curvature perturbation can be 
rewritten as 



7Z C =(1 - r)7Z$ + r7Z c 



=(l-r) 



n + _J_ ^(induced) 

1 — r 



=(l-r)7Z$ + r7Z a , (26) 

where we have redefined 71$ and 7Z a . After this redefinition, all the calculations presented in 
this section can be safely applied. 



4 Conclusions and discussions 

In this paper, we calculated the bispectrum of the curvature perturbation 7Z C produced during 
false vacuum inflation. Because 7Z C is intrinsically highly non-Gaussian, the corresponding 
bispectrum exhibits a novel shape. Using this result, we proposed the "old curvaton" scenario, 
where the Gaussian contribution of the total curvature perturbation consistent with the current 
observations is provided by the curvaton. The resulting shape of the bispectrum is a mixed 
type, i.e. the combination of the novel "false vacuum" and the well known local shapes. In 
addition to the novel bispectrum, in the old curvaton scenario we also find relatively large 
running of the spectral index. 

In the old curvaton scenario, both the inflaton and the curvaton play key roles in producing 
observable effects. Moreover, there are possibly thermal components to confine the inflaton at 
the top of the potential. These components set up interesting connections between our work 
and some other recent progress on inflation and non-Gaussianity: 

• Quasi-single field inflation: In the old curvaton scenario, the curvaton field slowly rolls 
down its potential. One could think of this scenario in an alternative way: to interchange 
the name of the curvaton and the inflaton, so that inflation is driven by the curvaton field a. 
Then our model describes a slow roll inflation direction, plus another direction along which 
the field is trapped at the minimum of the potential. In this viewpoint, our model looks 
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similar to quasi-single field inflation [12] . However, in quasi-single field inflation, the transition 
from the isocurvature perturbation to the curvature perturbation comes from direct coupling 
between the two directions. While in the old curvaton scenario, the transition comes from 
the non-conservation of curvature perturbation after the end of inflation, which is standard 
in the curvaton scenario. The origin of non-Gaussianity is also different. In quasi-single field 
inflation, non-Gaussianity originates from the interaction between the scalar fields. While in 
the old curvaton scenario, the non-Gaussianity originates from the non-linear mapping from 
the scalar field fluctuation to the curvature perturbation. 

• Thermal non-Gaussianity: In the old curvaton scenario, the effective potential fl3]) can be 
obtained in several ways. A simple way is that the inflaton is trapped at the origin for some 
e- folds by thermal radiations [22J. It is possible that they also have non-Gaussianities [23J, 
which could be transferred to the curvature perturbation by gravitational or other couplings. 
Here we note this possibility, and at the same time emphasize that this effect is under control 
in our model. Thus we can turn off these contributions, and our previous calculations give 
the leading non-Gaussianity. Nevertheless, it would be interesting to see whether there are 
parameter regions where these thermal effects can affect the running and the non-Gaussianity, 
as well as the power spectrum and the spectral index [24J. 

• Multi-stream inflation: At the end of false vacuum inflation, the potential (j3J) becomes 
unstable. Then, the field could roll towards either <p > or <fi < at different causal patches. As 
discussed in multi-stream inflation [25], these bifurcation does not necessarily lead to disasters, 
and can rather probably give rise to some interesting observational effects, such as features, 
non-Gaussianities, and the CMB asymmetries at the bifurcation scale. However, these effects 
do not change the calculation and conclusion in our present paper, because they are associated 
with perturbations on smaller scales. 

Finally, we would like to mention that the old curvaton scenario and the related works 
we discussed above all belong to the attempts assuming that the inflationary dynamics is 
more complicated than the simplest single field slow roll inflation. Recent progress in string 
theory, especially the string landscape [26] and related cosmology [27J, provides a theoretical 
playground, and also higher a priori expectation for the complicated inflationary dynamics. 
Hopefully, the old curvaton, together with a great number of other models and mechanisms, 
can serve as a building block for a realistic description of inflation in the string landscape. 
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A 7Z C in terms of (f> during false vacuum inflation 

In this appendix^, we provide a little more detailed argument on how we obtained @. We 
know from Ref. [TO] that the gauge invariant energy density perturbation on the comoving 
hypersurfaces pA is related to the gauge invariant intrinsic spatial curvature perturbation $ in 
the longitudinal gauge by a Poisson-like equation, 

pA = -2m| 1 (^) 2 V 2 $ . (27) 

To extract the functional behaviour of $, we can notice that from the energy density correlation 
function, 

D(x, x>) = AmUHvY <V^(*)V^(x')) ~ ^ , (28) 

so that $ has an exponential factor e _/ir . Since we are interested in the super-horizon sep- 
arations, the Laplacian operator primarily picks the terms with the least power of r in the 
denominator. Thus we can approximate 

V 2 $ « p 2 <5> . (29) 
Further, $ can be written in terms of the comoving curvature perturbation TZ C as 

$ = "32^ M Uc = -^rjP Uc ' (30) 

where A is a constant of 0(1), and (p + p) Tcn is the "renormalized" expectation value of p + p, 
given by 

(P + P)^ = A^- 2 ^t(p V ) 2 . (31) 

Therefore, combining fl29l) and fl30l . from (|27p we can write TZ C in terms of the density pertur- 
bation pA as 

f> W 2 . (32) 

We can move to the Fourier space more conveniently by noticing that for D(x, x') the most 
significant contributions come from the terms which do not contain any time derivative. This 
means, 



V 2 (pA) = V s 



£_ W) 2 

2 2a 



+ ^ + v(4>) 



+ 3#0V 2 



(V0) 
2a 2 



2 



+ V(<f>) 



(33) 



Further, as we are interested in the false vacuum inflation stage where m 2 s is positive, 1/a 2 
term in completely dominates the potential. Thus, after all, we have 

pA « ± [(V0) 2 + p 2 4> 2 ] • (34) 



7 Some parts of this section come from Ref. [17] . 
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Decomposing <p in terms of the Fourier mode 



d 3 k 



„ikx , 



we can easily find the Fourier component of pA as 



d 3 q 



2a? J (2tt 



,3^ 



ik-q| [v 2 + q ■ (q - k)] 



(35) 



(36) 



Substituting fl36|) into fl32|) . we can obtain 01]). 

Before we finish, we note that since the potential is quadratic, we can promote as a 
quantum harmonic oscillator and expand the Fourier mode <pk in terms of the annihilation and 
creation operators, and respectively, as 



[x 



d 3 k 

(>)3 



ik-x , 



■Xv) 



d 3 k 



ik-x 



a k y> k {r}) + aL k <ft(v) 



where and aj^. satisfy the canonical commutation relation 

[o k ,at] =(27r) 3 5( 3 )(k-q). 



(37) 



(38) 



The solution of the mode function (fik{v) satisfies the asymptotic behaviour (Q in the limit 
\v\ > !■ 
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